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The dynamics of star- forming regions — which 
mechanisms set the cluster formation efficiency? 
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Abstract 

The fraction of star formation that results in bound stellar clusters (cluster formation ef- 
ficiency or CFE) is a central quantity in many studies of star formation, star clusters and 
galaxies. Recent results suggest that contrary to popular assumption, the CFE is not (solely) 
set by gas expulsion, but is also influenced by the primordial environment, although its pre- 
cise behaviour remains unknown. Here it is discussed which mechanisms set the CFE, which 
recent advancements have been made to disentangle their contributions, and which studies 
are needed in the near future to achieve a quantitative understanding of the CFE. 



1 Introduction 



x 



Star clusters are popular tracers of the star formation process, both on stellar and galactic 
scales (e.g. |Smith et aT]|2007[ |Pflamm-Altenburg et aL] |2007t |Bastian et"aL]|2010|>. This 



approach assumes a certain fraction of star formation that produces bound stellar clusters, 
also dubbed the cluster formation efficiency (CFE). R has recently become clear that the 
CFE is not the product of two different 'modes' of dispersed and clustered star formation. 
Instead, star formation proceeds according to a continuous density spectrum, of which the 
high-density end yields bound star clusters and the low-density end consists of unbound 
associations ( Bressert et al.|20"l0 Gieles &: Portegies Zwart|2011 ). This suggests the existence 
of a (possibly varying) critical density that allows for the formation of bound structure. If 
the density spectrum of star formation would vary with environment, this would also imply 
that the CFE is environmentally dependent - dense star-forming regions should then yield a 



high CFE (see e.g. Adamo et al. 2011 for a discussion). 



2 Dynamics of star-forming regions 

It was originally thought that infant mortality (disruption by gas expulsion) would 
be the key mechanism behind a CFE lower than 100% (Lada & Lada 2003). The question 



thus arises how this new picture, in which stars do not form in clusters but in a hierarchical 
setting, connects to the concept of infant mortality, and also which other (environmental) 
processes may determine whether stellar structure ends up being bound. The following 
sections address the different mechanisms at play, and how Gaia can be used to estimate 
their relative contributions to the CFE. 



2 Infant mortality: internal disruption by gas expulsion 

In the classical picture of cluster formation, the removal of gas from the star-forming environ- 
ment by stellar winds and supernovae leads to a perturbation of the gravitational potential. 
If the gas fraction is large enough, this unbinds the cluster - the aforementioned 'infant mor- 
tality' (e.g. Bastian Sz Goodwin 2006 Goodwin &: Bastian||2006 ) . The indications for infant 
mortality are largely empirical. Surveys of young stellar clusters show that some 90% of all 



gas-embedded clusters does not survive the transition to the exposed phase (Lada & Lada 



2003). However, this need not imply that the decrease of the number of clusters is causally 



related to gas expulsion. 



As is shown by Bressert et al. (2010) and Gieles & Portegies Zwart (2011), not all stars 



form in clusters. At least a certain fraction of gas-embedded stellar structure is gravitationally 
unbound from the onset. While these unbound associations disperse on a crossing time, 
the bound part of the embedded structure could still be subject to infant mortality. The 
disruption of a cluster by gas expulsion requires a sufficiently massive central concentration 
of gas. The theoretical evidence for infant mortality is largely based on analytical estimates 
(e.g. Tutukov 1978; Hills 1980) or numerical experiments that obey such conditions - either 
by assuming a static g as potentiator dynam ical equilibrium between the gas and stars (e.g. 
Boily & Kroupa||2003MBaumgardt & Kroupa||2007|). 



It was recently found by Offner et al. (2009) that the velocity dispersion of stars in 



simulations of star formation are about a factor of five smaller than that of the gas. This 
result has recently been investigated further by Kruijssen et al. (2011a), where we analysed 



the dynamics of the stellar structure in the simulation by Bonnell et al. (2008). They identify 



subclusters using a minimum spanning tree (see Maschberger et al. 2010) and compute the 
virial ratios of the stellar component only, by ignoring the gravitational potential of the gas. 
This is equivalent to observing the system at the moment of instantaneous gas expulsion. The 
virial ratios are also corrected for binaries and higher-order multiple systems. The resulting 
distribution of virial ratios is shown in Fig. [T] for the ensemble population of subclusters from 
all snapshots. On average, the subclusters are very close to virial equilibrium when neglecting 
the gas, particularly after one free-fall time has passed. 

The virialised state of the subclusters is traced to the fact that they are generally gas- 
poor, with typical gas fractions below 10% within their half-mass radii. Since the simulation 



by Bonnell et al. ( 2008 ) does not include feedback by radiation or stellar winds, the low gas 



fractions must have a dynamical origin. About half of the gas depletion can be attributed 
to the accretion of gas onto the sink particles, while the other half is due to the accretion- 



1 Except for a normalisation of the gas potential that decreases with time when the gas is expelled. 
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Figure 1: Histogram of the virial ratios Q va of the subclusters from all snapshots of the 
simulation (solid line). The shaded histogram represents the set of subclusters from the 
last snapshot at t = 0.641 Myr (after approximately one free-fall time). The dashed line 
is a Gaussian fit to the data for all snapshots, with mean value Q v ir = 0.59 and standard 
deviation gq = 0.16. The vertical dotted lines again indicate the marginally gravitationally 
bound case (Q v ir < 1) and the virialised case (Q v ir = 0.5). 



induced shrinkage of the subclusters, which is the dynamical response of the subcluster to 



the mass increase of the sink particles (see Moeckel &; Clarke 2011). This indicates that the 
accretion of gas onto the sinks is sufficient to balance the overall gas inflow. The subclusters 
only become gas-rich (with gas fractions above 20%) beyond about three half-mass radii, and 
are thus embedded in an evacuated cocoon of gas. The lack of a central concentration of gas 
implies that the subclusters are relatively unperturbed by gas expulsion - a simple analytical 
estimate yields an ~ 8% expansion after the gas has been removed. 

The degree to which the subclusters are able to evacuate the surrounding gas dynam- 
ically depends on the free-fall time. For a certain density spectrum of star formation, a 
gas-poor state is more easily achieved in the high-density regions, which complete the largest 
number of free-fall times before the onset of gas expulsion. In low-density regions, the free- 
fall time is longer and the evacuation would be less efficient, implying that the disruptive 
effect of infant mortality could be more important. It is worth noting that the densities at 
which the effect of infant mortality is strongest thus coincide with the range where stars are 
predominantly formed in unbound associations. 



The cruel cradle effect: external disruption by the primor- 
dial environment 



The transition of young star clusters from the gas-embedded to the exposed phase is also 
accompanied by an environmental effect that decreases the number of clusters more strongly 
before and during gas expulsion than after it. Star clusters are tidally disrupted by passing 
giant molecular clouds (GMCs), particularly in dense environments, where the frequency and 
strength of the tidal shocks is high ( Gieles et al.||200"6 ). In recent work, it has been indicated 
that star-forming regions affect their offspring in this way: the GMCs in these regions are 
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Figure 2: Schematic representation of the relative importance of the mechanisms that set 
the cluster formation efficiency (CFE). Shown are the disruptive strengths of infant mortality 
and the cruel cradle effect as a function of ambient gas density. The slopes, normalisations, 
and the location of the cross-over point of the curves are all unknown. 



capable of efficiently disrupting the new-born stellar clusters (Elmegreen &; Hunter 2010 



Kruijssen et al. 2011b). In those cases where the ambient density is high enough, the tidal 



shocks would be capable of disrupting young clusters in a single encounter, independently of 
their mass (also see Gieles et al.|2006 ) . This form of enhanced disruption is most prevalent in 
star- forming regions because the mean disruption rate decreases as a cluster population ages. 



We identified the two responsible mechanisms for a decreasing disruption rate in Kruijssen 
et al. (2011b, see their Fig. 10): cluster migration, i.e. the motion of clusters away from 
dense star-forming regions into the field, and natural selection, i.e. the preferential survival 
of those clusters residing in less disruptive environments. Because the enhanced disruption 
terminates when the gas has been expelled, it influences the cluster population in a way that 
is very similar to infant mortality. However, rather than being an internal effect like infant 



mortality is, the primordial disruption by tidal shocks is external. In Kruijssen et al. (2011a), 
we therefore named it the cruel cradle effect. 

Contrary to infant mortality, which should peak in low-density regions, the cruel cradle 
effect decreases the CFE in high-density regions. This implies that in different environments, 
the CFE may well be determined by different physical mechanisms. In order to understand 



the relation between the CFE and the gas or star formation rate density (e.g. Goddard 
et al.|2010 ; Silva- Villa &; Larsen||2011 ; Adamo et al.||2011 ) , the relative contributions of these 
mechanisms need to be accounted for. The 'importance' of infant mortality and the cruel 
cradle effect is shown schematically in Fig. [2] as a function of the ambient gas density p am b- 
While the precise form of these relations is unknown, the evidence outlined above leads to 
the trends sketched in the diagram. Infant mortality becomes more effective as the ambient 
density decreases, with the curve flattening at very low densities to reflect the saturation that 
occurs when star formation is so dispersed that virtually all structure would be disrupted by 
gas expulsion. The efficiency of the cruel cradle effect increases with ambient density, with 
the curve potentially flattening at very high densities where (nearly) all newly formed clusters 
are immediately disrupted by the strong and numerous tidal shocks. The relation between 
the CFE and /9 a mb is then given by the product of the curves due to infant mortality and the 
cruel cradle effect. 
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Figure 3: Movie of two toy iV-body simulations to illustrate the kinematic difference between 
the disruption of a ~ 25 M (sub)cluster by infant mortality (left) and the cruel cradle effect 
(right). Infant mortality is modelled by removing a static background potential with the same 
properties as the (sub)cluster at t = 3 Myr, while the cruel cradle effect is characterised by 
the passage of a GMC with mass 1.25 x 10 4 M . The top panels show the projected spatial 
configuration of the stars in their centre of mass frame. In the top-right panel, the GMC 
is indicated with a diamond and a dotted circle, which denotes its Plummer radius. In the 
bottom panels, the directions of the velocity vectors of the stars are projected onto a sphere, 
with angles in degrees. Infant mortality retains the random orientation of the velocities, but 
the cruel cradle effect induces structure in velocity space. Click on the figure to start the 



movie. If your PDF reader does not support embedded movies, you can also access it here 



The relations in Fig.[2]apply to the fraction of star formation that takes place in initially 
bound (sub)clusters, meaning that star formation in unbound associations is not included. If 
the fraction of star formation taking place in unbound associations also exhibits a trend with 
ambient gas density, it should be multiplied with the result of Fig. [2] to obtain the actual 
CFE-p am b relation. In the coming years, combined theoretical and observational efforts 
should enable the quantification of the relative contributions of star formation in unbound 
associations, infant mortality, and the cruel cradle effect. 



Using Gaia to distinguish between the mechanisms that de- 
termine the cluster formation efficiency 



Because the early disruption of stellar structure due to infant mortality and the cruel cradle 
effect should occur on similar time scales, it is hard to distinguish between both mechanisms 
observationally. However, the kinematic properties of disrupted cluster remnants will depend 
on the process that destroyed them. When gas expulsion leads to the disruption of a young 
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cluster, the orientations of the velocity vectors of the stars remain random - after all, there 
are no external torques acting on the cluster. However, the cruel cradle effect has a different 
impact, since disruption by an external perturbation does lead to a preferential orientation 
of the stellar velocities. This is illustrated in Fig. [3j which shows a movie of two toy iV-body 
simulations in which a low-mass (sub)cluster is being disrupted. The striking difference in 
velocity space implies that it should be possible to distinguish between infant mortality and 
the cruel cradle effect. While Fig. [3] shows two highly idealised scenarios, it also indicates 
the key characteristics that can be used to observe in which environments infant mortality 
and the cruel cradle effect play a role. An additional advantage of considering the kinematics 
of young cluster remnants is that they need not be caught in the act of being disrupted. 
Depending on the local strength of the Galactic tidal field, the imprint of the disruption 
mechanism will remain visible for some time after the actual dispersal took place. 

The combination of radial velocities from spectroscopy with proper motion measure- 
ments and membership identification by Gaia will provide the 3D positions and space veloc- 
ities that are needed to reconstruct the processes shown in Fig. [3j Until these data become 
available, a theoretical effort should be made to consider the kinematics of star formation 
in unbound associations, infant mortality and the cruel cradle effect in more detail. This 
will expand the current qualitative understanding of the mechanisms that set the CFE to a 
quantitative census, which in turn will be essential for the predictive potential of star clusters 
as tracers of galaxy-scale star formation. 
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